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Comparison of seismic signatures of flares obtained by SOHO/MDI and 

GONG instruments 

S. Zharkov\ V.V.Zharkova^, S.A. Matthews^ 
ABSTRACT 



The first observations of seismic responses to solar flares were carried out using 
time-distance (TD) and holography techniques applied to SOHO/MDI Dopplergrams 
G^ . obtained from space and un-affected by terrestrial atmospheric disturbances. However, 

the ground-based network GONG is potentially a very valuable source of sunquake 
observations, especially in cases where space observations are unavailable. In this 
paper we present updated technique for pre-processing of GONG observations for 



p^. application of subjacent vantage holography. Using this method and TD diagrams 

Q . we investigate several sunquakes observed in association with M and X-class solar 

-1-^ . flares and compare the outcomes with those reported earlier using MDI data. In both 

ci . GONG and MDI datasets, for the first time, we also detect the TD ridge associated 
with the September 9, 2001 flare. Our results show reassuringly positive identification 

^ ! of sunquakes from GONG data that can provide further information about the physics 

OO ! of seismic processes associated with solar flares. 

r- 

Subject headings: Sun: photosphere; Sun:helioseismology; Sun:flares; Sun: oscilla- 
tions; Sun:data analysis 



yj 1. Introduction 

rS 

5^ ! Discovery of sunquakes associated with solar flares (i Kosovichev & Zharkovalll998l : lDonea et al 



19991) has opened a new era in the investigation of energy and momentum transport mechanisms 
from the upper atmosphere to the photosphere and beneath, uncovering structure of these spec- 
tacular events. Sunquakes, seen as circular or elliptical waves - ripples, propagating outward 
from impulsive hard X-ray (HXR) solar flares along the solar surface, appear 20-60 minutes after 
the flare onsets. The surface ripples are also associated with strong downward shocks preceding 
these ripples with close (1-4 minutes) temporal correlation with the start of HXR flares, indicating 
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that high energy particles pla y some role in initiation of sunquakes (|Zharkova & Zharkovl 120071 : 
Martinez-Oliveros et alJboOsk 



Even though every flare is expected to inject particle beams of one kind or another into a 
flaring atmosphere, inducing either shocks or magnetic impulses, not many of them have recorded 
measurable signatures of seismic activity. Initially only X-class flares were considered as candi- 
dates for producing sunquakes. The first flare detection used time-distance diagrams and reported 
well distinguished ripp les emanating from the center of the location of a hard X-ray source in the 
XI. 1 flare 9 July 1996 (iKosovichev & Zharkovalll998l : lDonea et al.lll999l) . It was followed by de- 
tection of quakes associated with two extremely powerful solar flares of class X17 and XIO which 
erupted in NOAA Active region 10486 on October 28 and 29, 2003. These two flar es, known as 
the Halloween 2003 flares, were extensively investigated in lOonea & Lindseyl OOOSb by applying 
subjacent vantage acoustic holography to SOHO/MDI Dopplergram data. The acoustic signatures 
were also shown to co-align with the hard X-ray signatures and GONG intensity observations re- 
vealed significant radiative emission with a sudden onset in the compact region encompassing the 
acoustic signature. 



Further analysis of the 29 October 2003 quake was presented in iLindsey & Doneal (|2008h . 
where the authors proposed a new method for correcting intensity data recorded by the Global 
Oscillation Network Gro up (GONG), that allowed comparison of acoustic k e rnels with white-light 
traces of the flare. Later Izharkova & Zharkovl (|2006r) : iKosovichevI (|2006[) : Izharkova & Zharkov 
(|2007h investigated the Halloween flare of 28 October 2003 using the time-distance diagram method 
applied to SOHO/MDI d ata and detected three dis tinct seismic sources that coincided with the 
holographic sources from lOonea & Lindseyl (l2005h. For the 29 O ctober 2003 flare there were 



no time-distance ridges found by the authors of IZharkova & Zharkovl (l2007b when they analyzed 



Kosovichev 



3). 



(20061). The first M-class flare 



MDI dopplergrams, however one was later reported in 
in which seismic signatures were detected by means of acoustic holography using Doppler veloc- 
ity data from SOH O/MDI instrument w as the flare that occurred in NOAA Active Region 9608 on 
September 9, 2001 doonea et al.ll2006al) . 



Later the list of acoustically active flares was significa ntly extended. iBesliu-Ionescu et al. 

(120051) reported another six such flares, adding another four in lOonea. A.-C. et al.l (|2006bb . In fact. 



during the period of observations with SOHO/MDI instrument there have been 17 flares showing 
signs of acoustic activity possibly related to sunquake^. This expanding list motivated researchers 
to look further and to explore the data from the ground-based GONG observatories, which offers 
better coverage of helioseismic data compared to SOHO/MDI. 

In the unusually quiet solar minimum between cycles 23 and 24 more attention was paid to 
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each new flare occurring on the Sun, one of which was the flare of 14 December 2006. The latter 
was observed only by the GONG instruments a nd revealed some notic eable seismic signatures in 
both time-distance ridges and egression powers ([Matthews et alJl201 lb . In order to validate these 
findings, a comparison is required of the signatures of sunquakes derived for both the time-distance 
and holographic techniques from GONG data with those from SOHO/MDI for a number flares 
with distinct seismic signatures. Such a comparison will allow us to understand the differences in 
appearances and to derive recommendations for the reliable detection of sunquakes from GONG 
data. 

In this study we consider three acoustically active flares that have the luxury of helioseis- 
mic observations available from both the GONG and the SOHO/MDI instrument. The available 
Dopplergrams (GONG and MDI) are used to analyze the ac oustic signatures of the flares b y us- 
ing both the time-distanc e diagram technique (TD method: iKosoyichev & Zharko vail 1998b and 
acoustic holography fe.g. lBraun & Lindseylll999l : lLindsey & Braunll2000b . The description of data 
and additional corrections for the technique applied to GONG data are presented in section [2l the 
results of the comparison are described in section [3] and conclusions with recommendations are 
drawn in section IH 



2. Description of data and techniques 



In this study we use three acoustically active flares with strong seismic signals detected by 
SOHO/MDI. The first flare is an M-class flare that occurred in NOAA Active Region 9608 on 
September 9, 2001. The solar quake associated with the flare was the first one dete cted for M-class 
flares and investigated by means of acoustic holography in iDonea et al.l (|2006ar ) using data from 
SOHO MDI. The other two are extremely powerful X-class flares that erupted in NOAA Active 
region 10486 o n October 28 and 29, 200 3, with associated sunquakes first detected using acoustic 
holography by iDonea & Lindseyl (|2005b . The Hallow een flares have also been in vestigated by 
applying time-distance method to SOHO/MDI data in Izharkova & Zharko vl (|2007b . where three 
distinct sources were detected for October 28 flare. H owever, after extensively an alyzing MDI 
velocity data for the 29 October 2003 flare, the authors of IZharkova & Zharko vl (|2007b did not find 
a time-distance ridge in any of time-distance diagrams computed around the flare location. 



2.1. Helioseismic techniques for quake detection 



In order to d etect and analyze the solar qu ake associated with a flare we use both time - 
distance analysis (JKosovichev & Zharko val 1 1998b and acoustic holography (JDoneaet al.lll999b 
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Time-distance analysis is applied to detect the circular ripples generated by the quake. This consists 
of rewriting the observed surface signal in polar coordinates relative to the source, i.e. v{r,9,t), 
and using azimuthal transformation 

p27T 

Vm{r.t)= / v{r,9,t)e-'^'d9, (1) 



to study the m = component for evidence of the propagating wave. Then if seen, the quake 
manifests itself as a time-distance ridge, thus providing estimates of the surface propagation speed 
and the time of excitation. In this work the GONG high-cadence velocity data were used in the 
time-distance analysis. 

Acoustic holog raphy is applied to calculate the egression power maps from observations. The 

holography method (IBraun & Lindseylll999l : lDonea et al.ll999l : lBraun & LindseylbOOdlLindsey & Braun 



2000 ) works by essentially "backtracking" the observed surface signal, ?/^(r, t), by using Green's 
function, G+(|r — r'\,t — f), which prescribes the acoustic wave propagation from a point source. 
This allow s us to reconstruct egress ion images showing the subsurface acoustic sources and sinks. 
Following lOonea & Lindseyl (l2005b . in temporal Fourier domain we have 



H+{r,iy)= / dh'G+{\r-r'\,u)^P{r',u), (2) 

where a, b define the holographic pupil and H^ (r , u) is the temporal Fourier transform of H^ (r , i) • 
Then 

//+(r,t)= / due'^^^'H^{r,u), (3) 



JAu 



whence the square amplitude of egression is called the egression power 



(4) 



Green's functions built using a geometrical-optics approach are used in this study. As flare acoustic 
signatures can be submerged by ambient noise fo r the relatively long period s over which the egres- 
sion power maps are integrated, again we follow lOonea & Lindseyl (120051) using egression-power 
'snapshots' to discriminate flare emission from the noise with pass-band integration in equation 
dS]) performed over positive frequencies only in order to reduce noise. Such a snapshot is simply a 
sample of the egression power within a time At = j^^ = 500s. The snapshots used in this work 
are taken from the egression power, P(r, t), at selected times, t. 



2.2. Observations and data reduction 



An M9.5-class flare occurred in NOAA Active Region 9608 around 20:40 UT on September 
9, 2001 at around 104° Carrington longitude and 26° latitude south. GOES soft X-ray flux reached 



-5- 



a peak at 20:46 UT, with the background emission remaining above the "C" level for most of the 
flare period. The flare of October 28, 2003, occurred in NOAA Active Region 10486 at 287° Car- 
rington longitude and 8° south latitude. It is classified as X17.2, one of the most powerful amongst 
the quake producing flares recorded. The GOES satellite detected increased X-ray flux starting at 
09:51 UT, reaching a maximum at 11:10 UT. On the following day in the same Active Region an 
XIO class flare occurred at around 270° Carrington longitude, 10° latitude in the southern hemi- 
sphere. The X-ray flux observed by GOES began to increase at 20:41 UT, reaching maximum at 
20:49 UT and ending at 21:01 UT. 



The he lioseismic observations analyzed in this study were obtained b y the GONG (IHarvey et al 



19881. 1 1996b ground-based observatories, and by the MDI instrument (jScherrer et al.l 1 1995b on 



board the SOHO spacecraft. Both GONG and MDI observe using the photospheric Ni I 6768 
A line. GONG observations are made with one minute cadence and normally include full-disk 
Dopplergrams, line-of-sight magnetograms and intensity images. In this work we use MDI full- 
disk Doppler images obtained with a one minute cadence. Velocity measurements in both cases are 
made from Doppler shifts of the Ni spectral line. MDI estimates the velocities from line instensi- 
ties (filtergrams) scanned in sever al locations across t he line, while GONG relies on a fast Fourier 



tachometric scans across the line ( Harvey et al.ll 19881) to derive the surface-velocity images based 



on a standard response of the line profile to Doppler motion caused by propagation of acoustic 
waves. 

For the October 28, 2003 flare we use two-hour-long full-disk velocity observations with one 
minute cadence from the SOHO/MDI instrument and GONG starting at 10:46 UT. The October 29, 
2003 and September 9, 2001 series commence at 20.00 UT on the corresponding dates. In addition, 
for the Halloween flares we use one-minute cadence intensity observations available from the 
GONG for the same period. Unfortunately, there are no such intensity observations available for 
the duration of the September 9 flare. Following the standard approach in local helioseismology, 
we extract datacubes centered on the region of interest from each full-disk series to remap the 
data onto the heliographic grid using Postel-projection and to remove a differential rotation at the 
Snodgrass rate. For the velocity data the series average full-disk velocity image is subtracted from 
each observations before the procedure, in order to remove the rotation gradient. Due to different 
resolution of the instruments, SOHO MDI data is remapped at 0.125 degrees per pixel resolution, 
while GONG datacubes are at 0.15 degrees per pixel. 

For acoustic holography we use Green's functions centered at 6mHz, so the datacubes are 
filtered in the frequency domain using a bandpass filter allowing the full signal in the 5-7 mHz 
band, with steep Gaussian roU-offs on each side. The pupil dimensions for each dataset for the 
selected flares are presented in Table [B 
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2.3. Additional corrections for the GONG data 



As GONG is a ground-based network, its data are affected by visibility conditions at the time 
of observation. Effects such as atmospheric smearing and local stochastic translation introduce 
spurious temporal variations in magnetized regions that can easily dominate over a genuine seismic 
signal. Another concern can be related to GONG usage of tachometric scanning of the Ni line, 
which, due to variations in atmospher ic conditions between the sta rt and end of the scan, is likely 
to be more affected (see, for example. iGrigor'ev & Kobanovll 1988b . 



The spurious Doppler shifts cited by Grigor'ev & Kobanovl (|l988l) are applicable to radiation 

incidence away from normal incidence (above 2°) passing through a Fabrey-Perot eta lon, which has 

effect ive path differences of 2.2 x 10^ wavelengths implemented for GONG (|Harvey & The GONG Instrument Tea 
1995b . The GONG optics ingeniously avoid this problem by directing the long optical path through 
glass and the shorter through air, the geometrical paths being the sam e to within about a micron 
(lTitle&Ram"seylll98d : iHarvey & The GONG Instrument Teamlll995b . 



In the presence of a strong magnetic field (e.g. sunspot u mbrae) Zeeman splitti ng of a mag 



netic line introduces spurious phase shifts in the measurements (|Rajaguru et al.ll2007 | ). One possi- 



ble reason for such effect is the reduced line intensities within a sunspot (|Toner & Labontdll993 



Brulsl I1993I : iNorton et al.l [2006) which cause noise such as due to variable atmospheric smear- 
ing to introduce spurious intensity o bservations from surrounding region into the desired pixel 
(lToner&Labontelll993l : lBrauJl997b . 



This, for instance, leads to significant differences in the computed acoustic power maps 
between the sunspot data from MDI and GONG, with the sp ace-based data gen erally showing 
suppression of the acoustic power over a sunspot region (e.g. lOizon et al.l 12009. and references 
therein), while the ground-based GONG data demonstrating a large power increase at the same 
location (see top row of Figure [I] for example) that is clearly noise-related due to the reasoning 
above. 

To correct the atm ospheric contribution in G ONG observations in the first instance we use 
the method developed in lLindsey & Doneal (|2008b where the intensity data are available, e.g. Hal- 
loween flares. The method works by measuring atmospheric seeing effects such as translation and 
smearing of GONG intensity observations in relation to a reference image and then removing their 
contribution from the intensity data. As both intensity and velocity data come from the same in- 
strument, we apply the parameters extracted from the intensity series to correct the line of sight 
velocity data. 

In addition, since the atmospheric noise affects mostly the measurements taken over magne- 
tized regions, we seek to minimize the contribution of such data to the quake- specific computation 
of egression power. First, we note that in the magnetized regions the atmospheric noise manifests 
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itself as spurious velocity fluctuations leading to a substantial increase of the observed acous- 
tic power. Such an increase in GONG data is assumed to be induced by the atmospheric seeing 
effects. Second, it is known that flares and associated with them solar quake sources are usually lo- 
cated over or near a sunspot. On the other hand, quake signatures such as ripples and time-distance 
ridges are normally seen in the surrounding non-magnetised region. This is, at least in part, due to 
the complex and less understood propagation of magneto-acoustic waves generated by the quake 
in the sunspot itself. Also, for our estimates of the egression power we use Green's function based 
on a non-magnetic model of the solar interior as it is intrinsic to the acoustic holography. In the 
view of equations ([SMl) it is then reasonable to minimize the contribution of the velocity data taken 
over magnetic regions (MRs) such as sunspot. 

For smaller sunspots this can be achieved by the choice of pupil, ensuring that the smallest 
radius is always selected outside of MR. When a sunspot is large, other methods will need to 
be considered such as weighting of the sunspot data in the velocity measurements. One possible 
option is to fully neglect such data, i.e. using zero as weights for sunspot velocities. This, of course, 
introduces artificial inhomogeneities in the computed egression power, but the qualitative strength 
of the quake source can still be evaluated by comparison with the egression power of surrounding 
plasma. In our experience, however, the best results are achieved by weighting all measurements by 
the inverse averaged acoustic power computed for the filtered velocity series. This is equivalent to 



normalizing the acoustic power of the filtered datacube, similar to the approach of iRajaguru et al 



(l2006h developed for time-distance helioseismology and phase-speed filtering. All of the GONG 
data used in the following sections are processed as described unless otherwise stated. 



Flare 


MDI pupil size 


GONG pupil size 


September 9, 2001 
October 28, 2003 
October 29, 2003 


15-60 Mm 
15-45 Mm 
15-45 Mm 


25-70 Mm 
25-95 Mm 
20-55 Mm 



Table 1: Holography pupil sizes for each dataset for the selected flares. 
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3. Comparison of the GONG and MDI results 

3.1. Comparison of acoustic holography results 

For illustration purposes, the results of calculations of the total egression power in 5-7 mHz 
range estimated from the MDI and GONG datacubes corresponding to observations of the Septem- 
ber 9, 2001 flare, are presented in the bottom row of Figure [B One can see that even after the 
corrections, the acoustic source suppression over the sunspot region is considerably weaker in the 
map computed from GONG data. This is generally the case for other observations we have consid- 
ered and is believed to occur because of the lower resolution and atmospheric noise contamination 
in the ground-based network's data. For these reasons, in order to reduce such contamination we 
consider larger pupil sizes when working with GONG data as shown in Table [B By choosing the 
larger lower limit on pupil dimensions we ensure the minimised contribution of the measurements 
taken over MR for egression power computation at the points near and around sunspots. 



3.1.1. Holography: September 9, 2001 

Egression power snapshots computed for the September 9 flare are presented in Figure [2l with 
the MDI data plotted in the left column and the GONG data on the right. Velocity images averaged 
over the series duration for both instruments (located at the top of the Figure [2l) demonstrate clearly 
the differences in the original data sets, which are du e to a further loss of resolution due to the 



atmospheric effects as described bv lLindsev & Doneal (l2008b . Our M DI egression measurements 



Doneaetal 



(|2006ah . Comparing these with the 



for this flare essentially duplicate the results of 

obtained GONG snapshots (Figure O, it is clear that even after the corrections, though many 

similarities are present, there is a significant variation between the two sets of images. 

The most obvious difference is the apparent absence of the region with low acoustic emission 
around the sunspot in the GONG produced data. Such an absence is clearly related to a much 
weaker signature of this phenomena in the GONG egression power seen in Figure [B This can 
be explained by the spurious atmospheric noise affecting GONG measurements over the regions 
with strong magnetic field. Nonetheless, the quake's signature is clearly present in the GONG 
measurements, with the locations of acoustic kernels agreeing very well for the two instruments. 
We note, however, the difference in acoustic kernel shape. This is, most likely, due to the reduced 
spatial resolution of GONG data suppressing higher-/ contribution to the egression. The possibility 
of atmospheric noise contamination is discussed later in Section IH 
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3,1.2, Holography: October 28, 2003 



Egression power snapshots computed for the October 28 flare are presented in Figure [3] with 
the reference GONG intensity image located in the top left corner, followed to the right by the 
MDI egression power snapshot taken at 1 1 :07 UT with the arrows pointi ng to the detected acousti c 



kernels. This image essentially duplicates the panel c) of Figure 3 in iDonea & Lindseyl (120051) . 
For a reference, the GONG snapshot for the same time is presented in the bottom row with and 
without arrows. Here, one can see the acoustic signatures at the same locations as in the MDI 
data. Again, we note the region with weaker lower acoustic emission as seen by GONG, which 
affects the visible contrast of the quake kernels relative to their surroundings. It is also clear that 
while the locations are the same, the shape and strength of each of the four kernels varies from 
one instrument to another. For example, source 1 (see Figure [3]) appears to be more prominent 
and extended in GONG measurements compared with MDI. Given the reservations about ground- 
based data outlined above with the fact that our correction procedure rather artificially modifies 
the oscillation amplitudes in GONG data, it is clear that MDI measurements are closer to the true 
picture of the event. Nonetheless, we reiterate that quake signatures are clearly visible in GONG 
measurements in the same locations as those detected by MDI. 



3,1,3, Holography: October 29, 2003 

Egression power snapshots computed for the October 29 flare are presented in Figured with 
the reference GONG intensity image located in the top left corner, followed to the right by the MDI 
egression power snapshot at 20:43 UT with the arrows pointing to the detected acoustic kernels . 



This image is essentially equivalent to the middle panel of Figure 6 in iDonea & Lindseyl (120051) . 
The corresponding GONG egression snapshot with the quake signature is plotted in the lower row 
on the left clearly present at approximately the same location as in the MDI plot. As an example, 
the egression power computed from the GONG velocity observation with the masked sunspot area 
is plotted in the bottom row to the right. 



3.2. Time-Distance diagrams 

3,2,1, 9 September 2001 flare 

The time-distance diagrams extracted from the MDI and GONG data are presented in Figure 
[51 The ridge representing the quake in the MDI image is relatively weak but can be clearly seen. 
As far as we know, this is the first time-distance ridge for a solar quake associated with an M-class 
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flare has been found. By comparing the MDI image with GONG time-distance diagram one can 
also detect in GONG a very similar disturbance located at the same part of the image. Although 
weaker and less defined than in MDI, nevertheless, the ridge is definitely present. Once again the 
relative weakness of the ridge can be explained by the fact that it is obscured by a significant noise 
contribution. As additional re-assurance, there is a near perfect coincidence between the MDI and 
GONG time-distance source locations. Also, as can be seen from Figure [2] where the location of 
time distance source is marked as plus sign on selected GONG plots, there is a good agreement 
between the egression acoustic kernels and time-distance source locations. 



3.2.2. Halloween flares 
We were not able to find any distinguishable time-distance ridges for the October 29 fl are in 



either MDI or GONG data, similar to the previous attempts fe.g. lZharkova & Zharkovll2007|) . The 
other GONG dataset for the October 28 flare, has a gap of about ten minutes between 11:30 and 
11:40 UT. However, we have attempted to build the time-distance diagram and the time-distance 
plot obtained from the interrupted GONG data. This is presented i n Figure [6l It shows ( at leas t 



a part of) a ridge, with the location corresponding to Source 1 in IZharkova & Zharkovl (|2007r) . 
reproduced in the top row of Figure [6l Once again, the location of the time distance ridge 
coincides with that of MDI and agrees with the egression measurements presented in section [3X2l 



4. General discussion and conclusions 

In this study we have compared egression power maps and time-distance diagrams derived 
from GONG and MDI data. SOHO MDI and GONG velocity datasets were used for three flares: 
M-class September 9, 2001 (Figures [MIS]), X-class October 28, 2003 (Figures [3] and [6]), X-class 
October 29, 2003 (Figure©. 

Reassuringly, the egression power map snapshots show seismic signatures common to both 
instruments for all flares. These signatures display an excellent agreement between the two instru- 
ments in terms of their time and location. We note, however, that even after the pre-processing, as 
outlined in Section [23l GONG egression measurements remain relatively noisy. This leads to im- 
portant differences from the MDI produced egression maps, which are only partially compensated 
by increasing the pupil sizes when working with GONG data. One such difference is the apparent 
variance in shape and strength of the detected acoustic kernels as seen by these two instruments. 
Another is the relative weakness of the suppression of acoustic sources below the sunspot. Since 
solar quakes are often located in the sunspot, this means that identifying such seismic signatures 
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in GONG egression measurements is a harder task due to its lower signal to noise ratio over the 
sunspot region. 

One such method of verification is the computation of the time-distance diagram. As we 
have demonstrated for September 9, 2001 and October 28, 2003 flares, such diagrams computed 
from GONG velocity data can present the additional evidence of the quake. Figures [5] and [6] clearly 
show that, in spite of being less sensitive, the GONG data can respond to the time-distance analysis 
producing noticeable ridges similar to those observed from the higher-quality MDI data. Results 
of the comparison with MDI time-distance measurements have again revealed an excellent spatial 
agreement between the two instruments in terms of the time-distance source location. Additionally, 
the fact that the locations of the sources observed with the GONG time-distance diagrams coincide 
with the acoustic kernels deduced from the GONG egression snapshots confirms that with these 
different techniques one observes the same events - seismic signatures induced by solar flares. 
Therefore, we conclude that the GONG data can respond to time-distance analysis. Obviously, 
due to the characteristic noise, not every quake can be expected to produce the time-distance ridge 
in GONG diagrams, but if a ridge is seen in the GONG ground-based data, one can expect that it 
will also be observed by using the higher-quality satellite MDI or HMI data. 

We believe, the results of this study show that quake detection based on helioseismic reduction 
of GONG observations is possible. However, as the data are subject to atmospheric smearing and 
other related instrumental effects, GONG observations clearly have less intrinsic sensitivity than 
the space-borne observations. A useful prospective object for further study might be the quanti- 
tative comparison of intrinsic sensitivities of ground- and space-based helioseismic observations 
under various seeing conditions. Nonetheless these results should allow us to add to the list of 
known sun-quakes by investigating the known flares in the Solar Cycle 23 using GONG data when 
MDI observations were not available. This will provide further information about the physics of 
seismic processes associated with solar flares. 
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Fig. 1. — September 9, 2001 M-class flare: comparison of unprocessed GONG and MDI integrated 
high frequency (5-7 mHz range) acoustic power maps (top row) and 5-7 mHz range egression 
power maps (bottom row) computed using MDI (pupil 15-60 Mm) and GONG (pupil 25 - 70 Mm) 
with correction for variable atmospheric smearing applied to GONG data. 



-13- 



MDI vel average 






GONC 


vcl avcraqc 












100- 










Lj 


# 


^ 


0- 


*« 


+ 




100- 


■M 


* 


< 



100 100 

Mm 



GONG egression power, 20:4, 




Fig. 2. — September 9, 2001 M-class flare: MDI data is in the left column, GONG is on the 
right, averaged MDI/GONG velocity i mage (top), followed by egression power snapshots. MDI 
egression plots reproduce the results in iDonea et al.l (|2006al) . Quake location is indicated by an 
arrow. The plus sign in two upper right frames indicates the source position assumed for the 
diagnostics specified by equation ([T]) presented in Figure El 
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Fig. 3. — October 28, 2003 X-class flare: Averaged GONG intensity imag e is top left, fol 
lowed by MDI egression snapshot to the right with arrows indicating sources in iDonea & Lindsey 
(120051) with #1, #3 and #4 corresponding respectively to time-distance sources #2, #1 and #3 in 
Zharkova & Zharkovl ([20071). At the bottom row: GONG egression p ower snapshot at 11:07 (le ft), 
and on the right is the same image with arrows indicating sources in IPonea & Lindseyl ([20051) as 
above. 
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Fig. 4. — October 29, 2003 X-class flare: (a) is GONG intensity image; (b) MDI egression power 
snapshot at 20:43; (c-d) GONG egression power snapshots around the quake time; (e) egression 
power computed from GONG velocity data using sunspot mask averaged over one hour. Location 
of the quake is indicated by an arrow; (f) egression power snapshot computed from masked GONG 
data used in (e) taken around the quake time. 
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Fig. 5. — September 9, 2001 flare: time-distance diagram computed from MDI velocity data (top 
row), and GONG dopplergram observations (bottom). 
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Fig. 6. — October 28, 2003 flare: (Top r6^w):time-distanc e diagram computed from MDI data. 
The plots are reproduced from lZharkova & Zharkovl (|2007h . (Bottom row): time-distance diagram 
computed from GONG data. along the j/-axis corresponds to 1 1:00 UT. 
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